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BeU-inequality violation and entanglement, measured by Wootters' concurrence and negativity, 
of two qubits initially in Werner or Werner-like states coupled to thermal reservoirs are analyzed 
within the master equation approach. It is shown how this simple decoherence process leads to 
generation of states manifesting the relativity of two-qubit entanglement measures. 
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I. INTRODUCTION 

Quantum nonlocality, responsible for violation of Bell- 
type inequalities 1, 2], and entanglement (inseparabil- 
ity) are the fundamental resources of modern quantum- 
information theory and still the most surprising features 
of quantum mechanics (see, e.g., Q)- It is therefore de- 
sirable to investigate the degrees of the BeU-inequality 
violation and of the entanglement of a quantum state 
not only in relation to efficiency of quantum-information 
processing, but also to understand better subtle aspects 
of the physical nature. 

It is well-known that pure states Q or a mixture of two 
Bell states violate Bell inequalities whenever they are en- 
tangled. So, one could naively think that the only mixed 
states that do not violate the Bell inequalities are sepa- 
rable states. However, Werner demonstrated the ex- 
istence of entangled states which do not violate any Bell- 
type inequality. The standard two-qubit Werner state is 
defined by 



p'^^O) = p\Y)(Y\ + 



1-p 



I® I 



(1) 



being a mixture for p E (0, 1) of the singlet state \Y) = 
(1 01) — \10))/^/2 and the separable maximally mixed 
state, given by / /, where / is the identity operator of 
a single qubit. Definition ^ is often generalized to in- 
clude mixtures of any maximally entangled state (MES) 
instead of the singlet state only. So, e.g^ one can analyze 
the Werner-like state defined by |^ S H 



p^)(0) = p\X){X\ 



l-p 



I® I, 



(2) 



as a convex combination of the Bell state \X) = (|00) + 
|ll))/\/2 and I ® I. The original Werner state (^1, con- 
trary to (O, is invariant if both qubits are subjected to 
same unitary transformation, say U ®U . Nevertheless, 
for a given p, both states (Q) and Q exhibit the same 
entanglement properties, thus 10) is also referred to as 
the Werner state 0, 0, la] ■ In addition to the fact that 



the Werner states can be entangled without violating any 
Bell inequality for some values of parameter p, they can 
still be used for quantum-information processing includ- 
ing teleportation 'ol Moreover, the Werner states, 
given by ^ and can be considered maximally en- 
tangled mixed states of two qubits 0| in the sense 
that their degree of entanglement cannot be increased by 
any unitary operations, and they have the maximum de- 
gree of entanglement for a given linear entropy (and vice 
versa) . 

We will study the effects of a lossy environment, mod- 
elled by thermal reservoirs, on the BeU-inequality viola- 
tion and on the entanglement of the initial Werner and 
Werner-like states in the quest for new states including 
those having different orderings induced by two entan- 
glement measures: concurrence C and negativity N (de- 
fined in Sect. IV). The relativity of the entanglement 
measures was first observed by Eisert and Plenio |12| |. 
They showed numerically, using Monte Carlo simulation, 
that the condition 



C{pi) < C{p2) ^ N{p^) < N{p2) 



(3) 



can be violated by some two-qubit mixed states, although 
it is satisfied if pi and p2 are the Werner states or pure 
states. Virmani and Plenio demonstrated that all 
good asymptotic entanglement measures are either iden- 
tical or impose different orderings of quantum states. 
The problem of relativity of the entanglement measures 
was also studied in Refs. [a, Eji HE Ea, Ii3 ■ Here, in par- 
ticular, we present analytical examples of different or- 
derings imposed by the concurrence and negativity for 
two-qubit states, which violate the Bell inequality to the 
same degree. 

The Letter is organized as follows. In Sect. II, we 
discuss a dissipative model and give a general solution 
for two decaying qubits initially in the Werner states. 
A comparative study of the BeU-inequality violation and 
entanglement of various decaying states are given is Sects. 
Ill and IV, respectively. A final comparison and conclu- 
sions are given in Sect. V. 



II. MODEL FOR LOSS MECHANISM 
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We analyze evolution of two initially correlated qubits 
subjected to dissipation modelled by their coupling to 



2 



thermal reservoirs (phonon baths) as described by the 
following Hamiltonian 



k.n 



k.n 



[5(/^)46w+h.c.: 



(4) 

which is the sum of Hamiltonians for the system, Hs^ 
reservoirs and the coupling between them, respectively. 
It is a prototype model, where qubits can be implemented 
in various ways, e.g., by single-cavity modes restricted in 
the Hilbert space spanned by the two lowest Fock states 
(see, e.g., Ilg)- In I0J, au is the annihilation operator 
for the kth [k = 1,2) qubit at the frequency w^; fol'^'* 
is the annihilation operator for the nth oscillator in the 
kth reservoir at the frequency f^i'^'' , and gi'^'* are the cou- 
pling constants of the reservoir oscillators. We assume no 
direct interaction between the qubits, thus the system 
Hamiltonian is simply given by Hs — /i^^^j^ Wfca^Ofe. 
The standard master equation for the model reads as 



d 



1 ^ 



k=l 



+ {nk + l)(2akpal - a\akP - pa\ak)} 



paka\) 
(5) 



where p is the reduced density operator for the qubits, 
7/. is the damping constant and is the mean num- 
ber of thermal photons of the fcth reservoir. The ex- 
act solution of © for arbitrary initial conditions and 
arbitrary-dimensional systems is well known. By confin- 
ing our analysis to the initial qubit states coupled to the 
quiet reservoirs {fii = ^2 = 0), the solution in the com- 
putational basis {|00), |01), |10), |11)} in the interaction 
picture can compactly be given as 



p{t) 



h 
h 
h 

V51.92P30 



ft 

hi 

V9192P21 

^/gi92P3i 



fi 

V51.92P12 

h2 



V9T92P03 

^/5^52Pl3 

5152/533 



(6) 

where the elements of the initial density matrix p(0) are 
denoted by p2t+k,2j+i = {i, k\p{0)\jj); gk = exp(-7fet) 
and 



fk = \/53-fe[(l - gk)p3,3-k + Pko] 

hk = 53-/c[(l - gk)p33 + Pkk)], 
h = 1 - (^1 + +5152^33)- 



(7) 



We will apply solution © to analyze the effect of dissi- 
pation on the Bell-inequality violation and entanglement 
of the initial Werner states. 



III. BELL-INEQUALITY VIOLATION 



We will study violation of the Bell inequality due to 
Clauser, Horne, Shimony and Holt (CHSH) In a 

special case of two qubits in an arbitrary mixed state p, 



one can apply an effective criterion for violating the Bell 
inequality: 



\Tt{pB, 



CHSH J 



< 2 



where Schsh is the Bell operator given by 



CHSH 



(b + b') 



(b-b') 



(8) 



(9) 



with its mean value maximized over unit vectors 
a, a', b, b' in 3?^. Moreover, cr is the vector of the Pauli 

spin matrices cri , 172 , cts , and scalar product a • cr stands 

3 

for ^ ijCTj. By noting that any p can be represented in 

the Hilbert-Schmidt basis as 

/ 3 
P = J i 1(E) I + r ■ a- ^ I + I (E)S ■ a- + tnm 

\ n,m— 1 

where r, s are vectors in 3?'^, Horodecki et al. [TqI |2C 
proved that the maximum possible average value of the 
Bell operator in the state p is given by 




max Tr (p^chsh) = 2 ^M{p) 

"cHSH 



(11) 



in terms of M(p) = maxj<fe {uj + Uk}, where Uj (j = 
1,2,3) are the eigenvalues of the real symmetric matrix 
Up — TJ Tp] Tp is the real matrix formed by the coeffi- 
cients tnm = Tr {pan ® Cm), and TJ is the transposition 
of Tp. Thus, the necessary and sufficient condition for 
violation of inequality © by the density matrix H1U|I and 
Bell operator (S_^for some choice of a, a',b, b' reads as 
M{p) > 1 (H¥(|. To quantify the degree of the Bell- 
inequality violation one can use M{p), 2 M(g)" (see, 
e-g-, Hill), or max{0, M{p) - 1 } (see, e.g., ,22j). But 
we propose to use the following quantity 



Bip) = v/max{0, M(p)-1} 



(12) 



which has a useful property that for any two-qubit pure 
state it is equal to the entanglement measures such as 
concurrence and negativity, defined in the next section. 
As M{p) < 2, it holds B{p) e (0, 1), where B{p) = 1 cor- 
responds to the maximal violation of inequality ((HJ and 
B{p) = for states admitting the local hidden variable 
model. The larger value of B{p) > the greater violation 
of the Bell inequality. Thus, B{p) can be used to quantify 
the degree of the Bell-inequality violation (BIV), and for 
short it will be referred to as the BIV degree. 

The BIV degrees for the initial Werner states Px\^) 
and Py\^) £^re the same and equal to 



S^P)(0) =max{0,2p2- l}i/2 



(13) 



implying that the states violate the Bell inequality iff 
l/\/2 < p < 1. By changing the parameter p into 
(f) = {1 — 3p)/2, 113|) goes into another well-known form 
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(see, e.g., [19|,|20|). The BIV degree of the maximally en- 
tangled states (p = 1) is also maximal and equal to one 
as shown in the left panels of Figs. 1, 2 and 4 at i = 0. 

By applying the solution JB)) one observes that the ini- 
tial Werner state p"^ (0) decays as follows 





(+) 



















(1 +P)5i52 

,(+) 



(14) 



where ^ {2~gi){2-g2)+pgig2 and 4 = 53-fc[2- 
(1 +p)gk]- By applying the Horodecki criterion one finds 
the eigenvalues of U to be ui_2 = P^ffi52 and ui^2 > 1*3 = 
[(1 — gi){l — (72) +Pffi52]^, which implies that the BIV 
degree evolves as 



B'i\t)^ms.x{0,2p^g,g2-1} 



1/2 



(15) 



Examples of decays of B^^\t) for p = 1 and 0.8 are pre- 
sented graphically by assuming that the damping con- 
stants are the same, for 71 =72 = 0.1 in Figs. l(a,b) 
and 2(a,b), or different, for 71 — 0.1 and 72 = 0, in Fig. 
4. If the qubits are initially in the Werner state py-*(0) 
then the evolution of the density matrix is described by 










4-) 

^2pV5i52 






'2^^5152 


-pgig2 and /il" 



(1 







P)9i92 
(16) 

53-fc[2- 



where ^ = (2-.gi)(2-52)-raiy2 ai^i^ 
(1 — p)gk]- Applying the same procedure as for the state 
px{t), one can find the following eigenvalues of U: ui^2 — 
p'^gig2 and ui,2 > "3 = [1 - (.9i + .92) + (1 - p)9i92]'^}- 
Although the third eigenvalue differs from that for px (t) , 
the BIV degrees Bx (t) and By (t) decay in the same man- 



B^^\t) = max{0, 2p'gig2 - 1}'^'' = B'^'it) (17) 

as shown, e.g., by solid curves in Figs. 2(a,b) and 4. 
Obviously, by changing the sign in the definitions of the 
Bell states \X) and \Y), one finds the same decay of the 
BIV degree. Thus, one could conjuncture that all MESs 
decay in the same way. We will show that this is not true 
by analyzing the following initial state 



|Z) = i(|00) + |01) + |10) 



|ll))^-i=(|0,+> + |l,-)) 



(18) 

where |±) = (|0) ± |1))/V2. State \Z) is another MES, 
which has the BIV degree (and also the concurrence and 
negativity) equal to one, and can be obtained from \X) 
by applying locally the Hadamard transformation to the 
second qubit. One can also define another Werner- like 
state as a mixture of the MES \Z) and the maximally 



mixed state given by / 

p'^i\0)=p\Z){Z\ 



I as follows (0 <p < 1): 
1-p 



-1(E) I 



(19) 



for which the BIV degree is given by (|13|l as for the stan- 
dard Werner state. For this and other reasons concern- 
ing the entanglement properties (see Sect. IV) being the 
same as for the state ^ , we shall simply refer to (|19f) as 
the Werner state despite the fact that ((T^ is not C/(8)C/ in- 
variant. The thermal reservoirs cause the decay of pl'^(O) 
as follows 



h' 

P9i^/92 
P\/gi92 



P91V92 P Vgi92 -P^/9192 
52(2-31) PV5152 -P^/9192 
Py^9i92 51(2-52) -^'5l^/52 
-P^j9i92 -pV3i52 -P9iV92 9i92 

(20) 

where h' = A + gig2 — 2{gi + 32)- The eigenvalues Uk are 
{p^5i52, ■^{v ± y^ti2 — Ap'^gfg2)} implying the following 
decay of the BIV degree 

B'i\t) = max{0,p2gig2 + l{v + ^v^ ~ Ap^gfgl) - l}i/2 

(21) 

where v = (1 - 51)^(1 - 32)^ + P^ 9192(91 + .92). For 
more explicit comparison of expressions for B^^'^y (0 and 
B^^\t), we find their short-time approximations up to 
linear terms in time for p > l/\/2 as follows: 



B^pit) = B'^'it) = max{0, Q - (71 + l2)^t + 0{t')}, 



(22) 



B<i^ (t) = max{0, Q - [5(71 + 72) - I71 " 72I] ^ + Oit')} 

(23) 

where Q — \/2p'^ — 1. By attentively comparing H17|l and 
(|^ . we find that it holds 



p^t 



B 



(p) 



it)^B'^\t)>B'i'it) (24) 

for any evolution times. In a special case of short times, 
this result immediately follows from H22() and H23() . So, for 
both the qubits coupled to the reservoir(s) (71,72 > 0), 

the evolution of B^^\t) > differs from that of B^^]y[t) 
as shown, e.g., by solid and broken curves in Fig. 2 (a,b). 
However, by assuming that only one of the qubits is cou- 
pled to the reservoir, say 71 ^ and 72 = 0, the BIV 
degree of p^^(i) decreases at the same rate as that of 

p'"^^ [t) and Py^ (t) and all the three the BIV degrees are 
given by 



b'-^^ (72 = 0, t) = B'^> (72 = 0,t) = B'i> (72 = 0, t) 

= max{0,2p2gi - l}i/2 (25) 

which is a special case of ((T7|) and ((?T)) . Clearly for 72 = 
0, H23I) goes over into H22() as expected. The decays of 
the BIV degrees for one of the damping constants equal 
to zero are presented graphically by solid curves in Fig. 
4. 



,ip)i 




FIG. 1: Decays of different measures for a given initial state in 
each sub-figure: The Bell-inequality violation degree B^\t), 
concurrence C^^^ (t) and negativity A''^''' (t) for qubits initially 
in the Werner states p^' (0) , where ip = X,Y, Z and p = 1 
(left panel) or p = 0.8 (right panel), coupled to the same 
reservoir with the damping constant 71 = 72 = 0.1. 

IV. ENTANGLEMENT 

To study the entanglement, we apply the concurrence 
and negativity being related to the entanglement of for- 
mation and entanglement cost, respectively. 

The entanglement of formation of a mixed state p is 
the minimum mean entanglement of an ensemble of pure 
states that represents p 2^: 

Ef{p)= mm Tp,E{\^P,) 0,^1) (26) 

where p = I]iP»IV'i)(^j| and i?(|V'i)(V'i|) is the entropy 
of entanglement of pure state {ipi) defined by the von 
Neumann entropy. As shown by Wootters [2J|, the en- 
tanglement of formation for two qubits in an arbitrary 
mixed state p can explicitly be given as 

Ep{p)^H[l[l + ^l-C{pr]) (27) 

where H{x) is the binary entropy and C{p) is the Woot- 
ters concurrence defined by 

4 

C(p) = max{0,2maxAj - V AJ (28) 

i — ^ 

i=l 



FIG. 2: Decays of a given measure for different states in each 
sub-figure. Notation and parameters are the same as in figure 
1. 




time 

FIG. 3: Negativities AN^/\t) = N'^^\t) ~ N^\t) {4< = 
X, Y, Z) corresponding to those presented in Fig. 2(f) and 
Table I for the same 71 = 72 = 0.1 and p = 0.8. 



p=1 p=0.8 




time time 

FIG. 4: Same as in figure 1 but for one qubit coupled to the 
reservoir with the damping constant 71 = 0.1, and the second 
qubit undamped (72 = 0). 
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where are the square roots of the eigenvalues of the 
matrix p{ay ay)p*{uy ® ay), where tXy is the Pauh 
spin matrix and the asterisk denotes complex conjuga- 
tion. Ep{p) and C{p) are monotonic functions of one 
another and both range from (for a separable state) 
to 1 (for a maximally entangled state), so that "one can 
take the concurrence as a measure of entanglement in its 
own right" pi |. 

The negativity is another measure of bipartite entan- 
glement being related to the Peres-Horodecki criterion 
iHlli and defined by [Hili^ 



N{p) = max{0, -2 ^ pj} 



(29) 



where p = pj^g is the density matrix of two subsystems 
(say, A and B with dA and cIb levels, respectively), and 
the sum is taken over the negative eigenvalues fij of the 
partial transpose p"^^ of p with respect to one of subsys- 
tems (say A) in the basis {|0), |1), I^a)}: 



cIa — I tie — 1 

E E 

i,j=0 kd=0 



(30) 



For two-qubit (d^ — dg — 2) states, the sum in 
can be skipped as p"^-^ has at most one negative eigen- 
value j30i |. The negativity, especially in low-dimensional 
systems (2 (g) 2 and 2 (g) 3), is a useful measure of entan- 
glement satisfying the standard conditions j2^ . The 
negativity (|29|l ranges from (for a separable state) to 1 
(for a MES) similarly to the concurrence and the BIV de- 
gree. It is worth noting that the logarithmic negativity, 
log2[iV(p) -I- 1], has a simple operational interpretation as 
a measure of the entanglement cost for the exact prepa- 
ration of a two-qubit state p under quantum operations 
preserving the positivity of the partial transpose (PPT) 

The concurrences and negativities for all the three ini- 
tial Werner states (0) (-0 = X, Y, Z) are the same and 
equal to 



<^(0) 



max{0,i(3p- 1)} (31) 



but different from their BIV degree B^^\o), given by 
(|13|l . For the Bell states, all the entanglement mea- 
sures and the BIV degree are equal to one. However, by 
decreasing parameter p, the BIV degree of the Werner 
states decreases faster than their entanglement. For ex- 
ample of p = 0.8, the initial values of the concurrences 
and negativities are equal to 0.7 while the BIV degree is 
0.529 as shown in the right panels of Figs. 1, 2, and 4. By 
comparing (|13|l and H31|l it is seen that the Werner states 
are entangled iffl/3<p<l. Thus, the Werner states 
for p e (l/3,l/-\/2) are entangled although admittinga 
local hidden model, i.e., satisfying the Bell inequality 

Qubits initially in the Werner state p^' (0) coupled to 
the thermal reservoirs exhibit dissipation described by 



(|14(l . We find with the help of the Wootters formula that 
the concurrence for p^^ (t) exhibits the following decay 

V9192 



C^f)(t) =max{0,. 



(2p - v/(2-gffi)(2-gff2)) } 

(32) 

with q = 1+p. While the negativity, according to Peres- 
Horodecki criterion applied for p^-* (t) , decays as 



N^Pit) = max{0,-[-gi-.g2 + (l+p)5i52 



V(.9i-.92)2 + 4p2gig2]}. (33) 



By assuming that both qubits are coupled to the same 
reservoir described the damping constant 7 = 71 =72, 
we observe that the concurrence and negativity are the 
same for all evolution times as described by 

d^\t) = N^^\t) = max{0, |[(1 +p).g - 2(1 -p)]}, (34) 

where g = exjp{—jt), as clearly depicted in Fig. l(a,b). 
In another special case, for the initial Bell states {p — 1), 
the solutions for the concurrence and negativity simplify 
to 

C^x\t) = V9^ (1 - \/(l-<?i)(l-<?2)) , 



N^xit) = 5132, 



(35) 



respectively. 

On the other hand, for qubits decaying from the initial 
Werner state Py\o), as described by (fT?^ . the concur- 
rence decays in time as 



1 



C^^'it) = max{0,-v/^(2p- 



xV(2-.gi)(2-52)-Wi52)} (36) 
while the negativity vanishes as follows 

N^''\t) = max{0,i(-2 + 5i+.g2-(l-p)gi.92 

+ v/(2-5i-52)2 + Vgi52)}. (37) 



Note that decays of all the entanglement measures, sim- 
ilarly to the BIV degree, are independent of the sign in 
definitions of \X) and In a special case for qubits 

coupled to the same reservoir (71 =72), and l|177|) 
simplify, respectively, to 

C^?) (t) = max{0, g (^p ~ ^ r{l - g) + \r^g^^ }, (38) 



max-IU, V(i - gf +p^g'^ ~ \rg^ - (1 - .9)} 



with r = 1 — p. In another special case for the initial Bell 
state {p = 1), the entanglement measures (|36|l and l|37|) 
reduce to 



(39) 
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respectively. It is worth mentioning that, in contrast to 
C^)(i) = N^^\t) for 71 = 72, the evolutions c'^\t) and 
NY\t) are different. Finally, let us briefly discuss the 

effect of dissipation of the initial Werner state p^^\o), 
described by (|20|l . on the entangled measures. The gen- 
eral analytical formulas for the p-dependent concurrence 
and negativity are quite lengthy thus are not presented 
here explicitly although were used for plotting the corre- 
sponding curves in Figs. l(e,f), 2(e,f) and 3. However, 
in a special case for the initial Bell-like state {p = 1) the 
decays of the entanglement measures are simply given by: 

C^z\t) = (l - iV(l-.9i)(l-52)) , (40) 

N^z\t) = y 9192(5 + 9192 - 2G) + 4(1 - G)2 + G - 1 

with G — (51 + .92)/2. By assuming the same damping 
constants 71 =72, e.g., the concurrence formula reduces 
to C'^P [t) — g{l + 9) ft- A graphical comparison of the 
decays of all the measures is shown in Fig. 2 and of the 
negativities in Fig. 3 and Table I for 71 = 72 = 0.1. 
Since the differences between the negativities n']^^ it) are 
not clear enough in Fig. 2(f), the curves were redrawn in 
Fig. 3 for the rescaled negativities AiV^''^ (<) = N^^^ (t) - 

N^\t) with = X,Y,Z. 

For clear analysis of the entanglement measures we will 
now focus on the special case for p = I. A comparison 
of Eqs. (123), ISni), and implies that the following 
inequalities are satisfied 

C'^\t)>c'i\t)>c'^\t) (41) 

for any evolution times. In the short time approximation, 
we find that the concurrences for p'"^^ {ip = X, Y, Z) decay 
up to linear terms in time as follows 

C^^it) = l-i(7i+72 + /^V^)t + 0(t2) (42) 

where /jf = 2, /y = 0, and = 1, which confirms 
the validity of (|41|l . Contrary to the concurrences, the 
negativities in the short-time limit decay as follows 

{t) = l- (71 + I2)t + \ [ll + 47172 + 72')^' + 0{t^) 

(43) 

where = 2, /y = 1, and — 5/4, which implies that 
at short evolution times the inequalities hold 

NP{t)>NP\t)>N':^\t). (44) 

A closer look at the inequalities if^ . (|^ and en- 
ables us to conclude that there exist states of two qubits 
(e.g., Px\t) and PyHt) at some short evolution time t) 
exhibiting the same BIV degree, B'^Pit) = B^it), but 
different entanglement measures in such a way that the 
concurrence C''-p{t) is smaller than Gy ''(t), while the neg- 
ativity N'^\t) is greater than NY\t). Obviously, the in- 
equalities 141|) for the concurrences correspond to those 



for the entanglement of formation, while the inequali- 
ties (|44f) for the negativities correspond to those for the 
PPT-entanglement cost. We stress that for longer times 
inequalities different from H44|l are satisfied for p < 1 
as presented in Table I. By analyzing this Table, other 
states differently ordered by the entanglement measures 
are readily recognized, including those at time t = t2, 
for which B^jf\t) = B^y^\t), N^x\t) = N^'^\t) but 

By comparing the series expansions (|42|) and (|43|l for 
71, 72 > 0, we can conclude that all the three negativities 
evolve in a more similar way (precisely they are the same 
up to linear terms in time) in comparison to more distinct 
evolutions of the corresponding concurrences. Also by 
analyzing (|22l) and (|43|l for p = 1, one observes that 

the negativities zi^) ^^"^ degrees _B^y(i) 

decay in the same manner up to linear terms in time 
for arbitrary values of 71 and 72. Other similarities of 
the decays of the entanglement and/or the BIV degree 
can be found for some special choices of the damping 
constants. In particular, H34|l is valid for 71 =72. In 
another special case of only one of the qubits coupled to 
the reservoir, e.g. 71 ^ and 72 = 0, the BIV degrees of 
all the three states decrease, as given by (|^ . but also 
the entanglement measures decrease at the same rate as 
expressed by their concurrence 

C"-^^ (72 = 0, i) = &p (72 = 0, i) = Cf (72 = 0, t) 
= max{0, i VsT (2p - V(l - P)[2 - (1 + p).9i]) } (45) 

and negativity 

N^f (72 = 0, <) = N^P'> (72 = 0, <) = N^'^ (72 = 0, t) 

= max{0, i (pgi + ^(1 - + 4pV - l)}- (46) 

By restricting to the case of the initial MESs {p = 1), the 
above equations reduce to {tp = X, Y, Z) : 

Cy)(72 = 0,t) = VsT, 
<^(72-0,i)=gi. (47) 

These conclusions are confirmed numerically on the ex- 
amples of 71 = 0.1, p = 0.8 and p = 1 as shown in Fig. 
4. 



V. CONCLUSIONS 

We have analyzed quantum-information properties of 
two decaying optical qubits prepared initially in Werner 
or Werner-like states and coupled to thermal reservoirs 
within the master equation approach. We have studied 
in detail a degree of violation of the Bell inequality due 
to Clauser, Horne, Shimony and Holt 0] by applying 
a parameter B related to the maximum possible mean 
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TABLE I; Comparison of the negativities for the three initial 
Werner states for p = 0.8, 71 = 72 = 0.1 as in Figs. 2(f) 
and 3. The characteristic evolution times are ti « 7.4745, 
t2 ~ 9.1613, ts ~ 9.5209. The corresponding inequalities for 
the BIV degrees are B^'(t) = -By'(i) > B^^\t), and for the 



concurrences are c'^\t) > C^'it) > C'^' {t) 



(p), 



time 



negativities 



t(p) 



> iV. 



t^ti 

t e (ti,t2) 

t = t2 

t e {t2M) 

t = ts 

t e (t3, 100) N^^>{t) > N'i>{t) > iV 



> ^. 
> 



X 



ifestations of the symmetry of our particular decoher- 
ence mechanism, the following properties of the decaying 
Werner states in relation to the Bell-inequality violation 
degree B and the entanglement measures C and N: If 
only one qubit is coupled to the thermal reservoir than 
those decays are independent of the initial Werner or 
Werner-like state for a given p. However, if both qubits 
are coupled to the reservoir (s) then the decays of C and 
N, and in some cases B, depend on the initial Werner 
or Werner-like state. By analyzing these decays, we have 
found states (say p and a) of two qubits exhibiting the 
same degree the Bell- inequality violation, Bp = B„, but 
different entanglement measures in such a way that the 
concurrence Cp is smaller than Ca, while the negativity 
Np is greater than Na- We have also found other states 
p and (7, for which either (i) Bp = B„, Cp < and 
Np< N„, or (ii) Bp ^ B„, Np ^ N„ and Cp < C^. Thus, 
the analysis of the decaying Werner states shows clearly 
the relativity of two-qubit entanglement measures. 



value of the Bell operator in a given state according to 
the Horodecki criterion On the other hand, the 

degree of the entanglement was expressed by the Woot- 
ters concurrence C as a measure of the entangle- 
ment of formation, and by the negativity N based on 
the Peres-Horodecki criterion '25*, z6\ and related to the 
PPT-entanglement cost ^J. We have observed, as man- 
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